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Recent observations^ 2, 3 have revealed an unexpectedly high 
binary fraction among the Trans-Neptunian Objects (TNOs) 
that populate the Kuiper Belt. The TNO binaries are strikingly 
different from asteroid binaries in four respects^: their frequency 
is an order of magnitude larger, the mass ratio of their compo¬ 
nents is closer to unity, and their orbits are wider and highly 
eccentric. Two explanations have been proposed for their for¬ 
mation, one assuming large numbers of massive bodies^, and one 
assuming large numbers of light bodies^. We argue that both as¬ 
sumptions are unwarranted, and we show how TNO binaries can 
be produced from a modest number of intermediate-mass bodies 
of the type predicted by the gravitational instability theory for 
the formation of planetesimals^. We start with a TNO binary 
population similar to the asteroid binary population, but sub¬ 
sequently modified by three-body exchange reactions, a process 
that is far more efficient in the Kuiper belt, because of the much 
smaller tidal perturbations by the Sun. Our mechanism can nat¬ 
urally account for all four characteristics that distinguish TNO 
binaries from main-belt asteroid binaries. 

The TNO binary 1998WW31 has^ a mass ratio m 2 /mi ~ 0.7, eccentricity e ~ 0.8, 
semimajor axis a ~ 2 x 10^ km, and inferred radii ri ~ l.lr 2 ~ 10^ km, hence 
a/ri > 10^, in stark contrast to main belt asteroid binaries'^, where m 2 /mi -C 1, 
0 , and a/ri < 10. 

Asteroid binaries are probably formed by collisions®, as in the leading scenario for 
the formation of the Moon^’ xhe observed characteristics, 1712 /mi -C 1, e ~ 0, and 
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a/ri < 10, are all natural consequences of this scenario^i. For a different scenario 
for 1998WW31, we can look at dynamical binary formation in star clusters, where 
there are three channels: 1) tidal capture^^j 2 ) three-body binary formation^^j ^^d 3) 

exchange reactions^3_ 

Channel 1 is analogous to the standard scenario for asteroid binary formation. It 
will indeed occur: each TNO has grown through accretion, and much of this accretion 
has happened through collisions with an object comparable in mass to the growing 
TNO itselfi4,15. 

Channel 2 would require a near-simultaneous encounter of three massive objects 
with low enough velocities to allow an appreciable chance to leave two of the objects 
bound. For this to work, the random velocities of the most massive objects should 
be signihcantly lower than their Hill velocities. Under such conditions, this channel 
could play a role, as pointed out by Goldreich et al.^, who assumed that there are 
~ 10® 100 km-sized object embedded in a sea of small (< 1km) objects. This as¬ 
sumption, however, is at odds with Goldreich and Ward’s theory for the formation 
of planetesimals® through gravitational instability, and it is hard to see how objects 
in the Kuiper Belt could form from non-gravitational coagulation, because the time 
scales are far too longi6. In contrast, the gravitational instability theory predicts the 
size of the initial bodies to be 10 — 100km. Starting with these larger bodies would 
make channel 2 ineffective, because the velocity dispersion would be higher than the 
Hill velocityi4. i7. 

Recently, Weidenschilling^ proposed a variation on the idea of using interactions 
between three unbound bodies in order to create a binary. He studied how a third 
massive body could capture the merger remnant from a collision of two massive 
bodies if the third body were near enough during the time of the collision. This 
mechanism seems unlikely to work, however, since it requires a number density of 
massive objects about two orders of magnitude higher than the value consistent with 
present observations^. 

Goldreich et al.^ have proposed another mechanism, based on the dynamical fric¬ 
tion from a sea of smaller bodies that can turn a hyperbolic encounter between two 
massive bodies into a bound orbit under favorable conditions. Effectively, this mech¬ 
anism makes use of a superposition of three-body encounters, since each light body 
interacts independently with the two heavier ones, and in that sense it is another 
variant on channel 2. As we mentioned above, the gravitational instability theory for 
the formation of planetesimals® would exclude the existence of such a sea of small 
objects, and since the alternative theory of nongravitational agglomeration does not 
seem to work, we will explore the consequences of dropping channel 2. 

Ghannel 3 can operate on the binaries formed through channel 1, so we should 
check whether channel 1 and 3 together produce the right binaries in the right num¬ 
bers. 

Starting with the hrst task, consider a relatively massive TNO primary in a binary 
orbit with a much less massive secondary. If the binary encounters a particle with a 
mass m that is comparable to the mass of the primary component (mi ~ m ^ ^ 2 ), 
the most likely result is an exchange reaction, in which the incoming object replaces 
the original secondary's. Figure 1 shows an example of such a reaction. 

The binding energy of the binary will not change much during the exchange, 
hence mim 2 /ao ~ rriim/a where a is the new semimajor axis after the exchange. 




Figure 1: An example of a binary-single-body exchange interaction, in the ‘(massive, 
light) meets massive’ category discussed in this paper. Bodies 1 and 2 have masses 
rrii = 1 and m 2 = 0 . 1 , respectively, forming a binary with an initially circular orbit. 
Body 3, with mass m 3 = 1, encounters the binary on an initially parabolic orbit. In 
panel (a), the whole scattering process is shown. Panel (b) shows the complex central 
interaction in more detail, while panels (c) and (d) show the orbits of the initial and 
hnal binary, respectively. Note that the hnal binary orbit is highly eccentric and 
much wider than the initial circular binary orbit. 
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Table 1: Cross sections a for various configuration-changing channels in binary-single¬ 
body scattering. The gravitational focusing factor is scaled out in order to obtain 
hnite values in the parabolic limit, where v is the initial relative velocity between 
binary and single body at inhnity. We use units in which G = mi = m 3 = a = 1, 
where G is the gravitational constant, mi and m 3 are the masses of the heaviest body 
in the binary and the single body, respectively, and a is the initial semi-major axis of 
the binary. The mass of the lighter body in the binary is m 2 = 0.05. The radii are 
ri = r 3 = 0.05 and r 2 = ri(m 2 /mi)^/^ ~ 0.01842. The scattering processes are coded 
as follows: {x, y) indicates a binary in the hnal state with components x and y, while 
a p + q indicates the product of a merger between bodies p and q. A single body 
2 ; in the hnal state is indicated by {,),z. The physical meaning of the six channels 
is as follows: (a) an exchange reaction resulting in a massive-massive binary; (b) 
an exchange reaction resulting in a massive-light binary; (c) a merger resulting in a 
massive-massive binary; (d) a merger resulting in a twice-as-massive-light binary; (e) 
a merger resulting in a massive-massive binary; (f) no binary is left, after three-body 
merging or two-body merging followed by escape. 


This implies a/ao ~ m/m 2 ^ 1. Under the impulse approximation, the interaction 
happens in a space small compared to the distance oq to the primary. Conservation 
of specihc angular momentum of the system gives m 2 ao(l — Cq) ~ ma(l — e) which 
gives 1 — e ~ m 2 /m <C 1. 

We have run a series of scattering experiments to obtain the relevant cross sections, 
for an initial binary with mass ratio of 20 : 1 and semimajor axis Oq = 20 ri, where ri 
is the radius of the primary. These values are typical for main-belt binary asteroids, 
with m 2 /mi < 0.1, and separations 5 — 40 times the radius of the primary. We choose 
parabolic relative orbits for the single body approaching the binary, with periastron 
distances uniformly distributed between 0 and 20ao. We only followed the system as 
long as all three bodies stayed within their Hill radius, lOOOao. 

Table 1 gives cross sections for processes in which initial binary membership is 
altered. Channels (a), (c) and (e) result in binaries with two massive components, 
and together comprise about 80% of the total cross section. We checked these results 
through a comparison with the starlab three-body scattering package^^. 

In hgure 2 the distribution for the semi-major axis is strongly peaked at a = 20, in 
good agreement with the simple argument presented above. Similarly, the eccentricity 
peaks at 0.95, as expected. We assumed ri = 75km, the estimated radius of the 
primary of 1998WW31. In hgure 3, the orbital elements of 1998WW31 are consistent 
with the binary having formed through the processes modeled here. 

We now confront our second task: to check whether the exchange channel is 
efficient enough to produce the observed binaries. Starting with TNOs of inter¬ 
mediate mass, as predicted by Goldreich and Ward’s theory for the formation of 
planetesimals®, the heaviest TNOs will accrete mass primarily through collisions with 
TNOs of comparable mass^"^’ Many of these collisions are of the ‘giant impact’ type 




a[km] 

0 2x10"^ 4x10"^ 6x10"^ 



e 


Figure 2: Normalized differential cross sections for the formation of a ‘massive- 
massive’ binary, under the conditions specihed in the text (channels a, c and e in 
table 1), with respect to the semi-major axis a (top panel), and eccentricity e (bottom 
panel) of the hnal binary. The initially circular binary has a = 1 in the dimensionless 
units used for da/da, while the physical units are given for reference at the top of the 
figure. The hlled points are the total values for the differential cross sections, while 
the open circles are the contributions from the merger channels (c and e in table 1). 
Note the double-peaked structure in the top panel: the sharp peak toward a ~ 20 
arises from non-resonant exchanges, where the hnal binary has an energy comparable 
to that of the initial binary; the broad peak around a ~ 10 arises from resonant 
exchanges, where the memory of the initial binary is wiped out, leading on average 
to more strongly hyperbolic escape in which a harder binary is formed. 
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Figure 3: Orbital properties of ‘massive-massive’ binaries formed in our scattering 
experiments: a and e have the same meaning and units as in fig. 2. Contributions 
from exchange reactions, channel (a) in table 1, are limited by energy conservation to 
a < 20, and give rise to the horizontal rim in the middle of the hgure. Contributions 
involving mergers, channels (c) and (e) in table 1, can lead to a values all the way 
to the Hill radius a ~ 10^, but are limited by angular momentum conservation to 
increasingly high e for increasing a. The star symbol shows the observed orbit for 
1998WW31. Boxes around the star indicate the observational 1- and 2-a error bars. 

that form a tight circular strongly unequal-mass binary (channel 1). Let us estimate 
what fraction of encounters between comparable-mass TNOs will give rise to ‘giant 
impact’ type binaries, and how long such binaries survive on average before they are 
destroyed again. 

We assume that one in three collisions between comparable TNOs gives rise to a 
binary. When no binary is produced, we have to wait for a typical time T until another 
major collision occurs. When a binary is formed, gravitational focusing implies a 
cross section for three body interactions of order oq. Therefore, our newly-formed 
binary will undergo an exchange reaction on a time scale (r/ao)T <C T, leading to a 
signihcant increase in a. Strong three-body interactions will subsequently occur on 
a much shorter time scale {r/a)T T. As a result, the semimajor axis will shrink 
systematically, while the ‘thermal’ distribution /(e) = 2e favors high eccentricity^^. 

When the orbit becomes small enough, r/a ~ 0.03, the chance for collisions in 
resonant encounters becomes signihcant^o. Let us assume that an exchange reaction 
turns a ‘giant impact’ binary into a binary with a semi-major axis of a ~ 300r. Each 
subsequent strong encounter will on average decrease a by a factor^i ~ 1.2. After a 
dozen encounters, a ~ 30r and collision is likely to occur. The time scale for each 
encounter to occur is ~ {r/a)T. The waiting time for the last encounter in this 
series to occur is (1/30)T, while each previous waiting time was less by a factor 1.2. 
Summing this series, we get a total waiting time of (T/30)/(l —(1/1.2)) = 0.2T before 



a collision between two or three massive TNOs. If all three collide, we are back where 
we started, and the resulting system may be a single body (with an assumed chance 
of 2/3) or a strongly unequal-mass binary (chance 1/3). If two of the bodies collide, 
the third one may remain in orbit, or it may escape. In the latter case, we again are 
back where we started. In the former case, we still have an equal-mass and likely 
highly eccentric wide binary. 

Under these assumptions, in 1/3 of the cases, we wind up with an equal-mass 
TNO binary with the observed properties for a period ~ 0.2T, compared to a 2/3 
chance to wind up with a single TNO for a period ~ T. This allows us to derive the 
rate equation for the formation and destruction of the binaries. If we denote by Ns 
and Nb the number of single bodies and the number of binaries, respectively, we have 


dNe 

dt 

dNs 

dt 


-Ns - — -Nb 
3 0.2 3 

dNB 

dt 


if we measure time in unit of T. So for the stationary state we have dNB/dt = 
dNs/dt = 0, and Nb = 0 . 2 iV 5/2 = O.lNs- Therefore, the binary fraction is ~ 10%. 
When accretion in the Kuiper belt region diminished, the number of single and binary 
objects was frozen, with a ratio similar to this steady-state value. 

While our arguments are only approximate, it is clear that after cessation of the 
accretion stage at least several percent or more of the TNOs were accidentally left in 
such a binary phase. The fact that more than 1% of the known TNOs are found to 
be in wide roughly equal-mass binaries is thus a natural consequence of any accretion 
model independent of the assumed parameters for the density and velocity dispersion 
of the protoplanetary disk or the duration of the accretion phase. As a corollary, 
we predict that future discoveries of TNO binaries will similarly show roughly equal 
masses, large separations, and high eccentricities. 

We conclude that we have found a robust and in fact unavoidable way to produce 
the type of TNO binaries that have been found, as long as we start from the plausible 
assumption that TNOs were formed through gravitational instabilities®. 


1. Burnes, J. A., Science, 297, 942-943 (2002). 

2. Veillet C., et al.. Nature, 416, 711-713 (2002). 

3. Margot, J-L., Nature, 416, 694-695, (2002). 

4. Weidenschilling, S. J., Icarus, 160, 212-215 (2002). 

5. Goldreich, P., Lithwick, Y., and Sari, R., Nature, 420, 643-646 (2002). 

6. Goldreich, P. and Ward, W. R., ApJ, 183, 1051 (1973). 

7. Margot, J-L., Nolan, M. C., Benner, L. A. M., Ostro, S. J., Jurgens, R. F., Giorgini, 
J. D., Stade, M. A., and Campbell, D. B., Science, 296, 1445-1447, (2002). 



8. Merline, W. J. et al. in Asteroids III (eds Bottke, W. F., Cellino, A., Paolicchi, P. 
Binzel, R. P.) (in the press). 

9. Hartmann, W. K. and Davis,D. R., Icarus, 24, 504-515 (1975). 

10. Kokubo, E., Ida, S., Makino, J. Icarus, 148, 419-436, (2000). 

11. Durda, D. D., Bottke, W. F., Asphaug, E., and Richardson, D. C., Bull. Am. Astron. 
Soc., 33, 1134 (2001). 

12. Fabian, A.C., Pringle, J.E. Rees, M.J., MNRAS 172, 15P (1975). 

13. Heggie, D.C., MNRAS, 173, 729 (1975). 

14. Kokubo, E. and Ida, S., Icarus, 131, 171 (1998). 

15. Makino, J., Fukushige, T., Funato, Y., Kokubo, E., New Astronomy, 3, 411-416 
(1998) 

16. Wetherill, G.W., Annu. Rev. Earth Planet. Sci. 18, 205 (1990). 

17. Kenyon, S. J.and Luu, J. X., Astronomical Journal, 115, 2125 (1998). 

18. Spitzer, L., "Dynamical evolution of globular clusters”, Princeton, NJ, Princeton 
University Press, (1987). 

19. McMillan, S. Hut, P., Astrophysical Journal, 467, 348-358 (1996). 

20. Hut, P. and Inagaki, S., Astrophysical Journal, 298, 502-520 (1985). 

21. Heggie, D.C. Hut, P. The Gravitational Million-Body Problem [Cambridge Univ. 
Pr.], Ch. 23 (2003). 

Acknowledgements We acknowledge helpful comments on our manuscript by 
Peter Goldreich and Roman Rafikov. 



